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Abstract 

We study the dynamics of a charged particle in a planar magnetic field which consists of 
n > 2 disjoint localized peaks. We show that, under mild geometric conditions, this system 
is semi-conjugated to the full shift on n symbols and, hence, carries positive topological 
entropy. 

1 Introduction and results 

Consider a magnetic field in R 3 whose field lines are perpendicular to the plane R 2 x {0} = R 2 . 
Then the motion of a particle of unit mass and unit charge in that plane is modelled by 
Newton's Second Law 

q = B(q)Jq (1) 

where B : R 2 — > R describes the field strength and the term on the right hand side is 
the Lorentz force corresponding to the magnetic field, with J being the symplectic matrix 
( - i o ) ■ The differential equation (1) can be written as the Hamiltonian system generated by 
the Hamiltonian H : T*R 2 ->• R,H(q,p) = ±\\p\\ 2 on (T*R 2 ,w) with the twisted symplectic 
form 

lu = loq + B(q)dqi A dq 2 

were ujq = d\ = dpi A dqi + dp% A dq 2 stands for the standard symplectic form on T*R 2 . 

In this paper, we study the dynamics of a particle when the magnetic consists of n > 2 
localized peaks, i.e., when the support of B consists of n connected components suppSfc. 
Assuming that each component is a disc where the magnetic field is rotationally symmetric, 
we can show that this dynamical systems exhibits chaotic behavior in the sense that it 
possesses an invariant set on which the Poincare map induced by its flow is semi-conjugated 
to the full shift in n symbols. This implies that there are solutions visiting the different 
components suppBfc in any prescribed order; see Fig. 1. 

Moreover, we can conclude that for all sufficiently small energies E > our system has 
topological entropy h top > cVE > 0. 

There is a vast literature on similar results for the motion of a particle in a multi-bump 
potential where the Hamiltonian is given by H(q,p) = ^||p|| 2 + V(q), and uj = wo- Classical 
scattering by potentials is treated in [DC]. References on chaotic scattering can be found, 
e.g., in [Sm, RR]. Scattering by obstacles is discussed in Sect. 5 of [Ga]. Many of these works 



1 



supp Bi 
/ ^=\~^ v supp B 2 




suppS 3 



Figure 1: Symbolic Dynamics: Trajectory realising the prescribed order .., 1, 2, 3, 3, 2, .. 

(including the case [Kn2] of molecular scattering) contain an analysis of the bounded orbits 
using symbolic dynamics. In the abovementioned cases the base for such an analysis is that 
scattering by a single center is described by a nontrivial topological index [Knl, KK]. 
Recent progress on motion in magnetic fields includes the following results. 

• It is known that the horocycle flow on the upper half plane can be interpreted as a 
magnetic geodesic flow. Dividing out co-compact subgroups of PSL(2,R), this gives 
rise to magnetic geodesic flows on compact surfaces without periodic orbits. However 
both the metric and the magnetic field where shown in [Pa] to be rigid in this case. 

• This was used in [Gi] to provide counterexamples to the Hamiltonian Scifcrt conjecture 
for M. 2d , d > 3 (that is, constructing a smooth proper Hamiltonian exhibiting a regular 
level set without periodic orbits). 

• In [Gr] bounds for entropies of the magnetic geodesic flow on negatively curved man- 
ifolds were obtained, and topological entropy was shown to be strictly larger than the 
Liouville metric entropy except in the constant curvature and zero field case. 

For magnetic geodesic flows on general closed surfaces criteria for positive topological 
entropy were derived in [Mi]. Relationships between the growth rates of magnetic 
geodesies and the topological entropy can be found in [Ni] . 

• In [PS] the magnetic geodesic flow was compared with the geodesic flow for energies 
above the Mane critical value. Mane's action potential was found to coincide with the 
Riemannian length. 

Surprisingly, the study of the dynamics of multi-bump magnetic fields has, to the best 
of our knowledge, not been undertaken yet. Note that, in contrast to the case of a singular 
potential, the magnetic Hamiltonian flow is complete. Indeed, since the energy E = H is 
constant along the trajectories, we have 

\\p(t,x )\\=V2E 

and 

||g(t,a:o)|| = ||g(0,a;o) + /oP(*,a:o)|| < ||?(0, x )\\ + \t\ V2E 
for all t £ K. Moreover, the energy surfaces := H^ 1 {E) with E > are all diffeomorphic 
to IR 2 x S 1 , so the embedding of symbolic dynamics will not be based on the nontrivial 
topology of the energy surface but needs to be established by other methods. 
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The paper is organized as follows. In Section 2, we will first study the dynamics generated 
by one component, i.e., the dynamics in a rotationally symmetric magnetic field. The main 
results will be then stated and proven in Section 3. Finally, in Section 4, we make some 
further remarks and point out open problems we will attack in the future. 

Acknowledgements: The second author acknowledges financial support by the German 
National Academic Foundation. 

2 Dynamics in a rotationally symmetric magnetic field 

First of all, it is necessary to understand the dynamics generated by one of the components. 
Therefore we consider the case of a compactly supported, rotationally symmetric magnetic 
field B : R 2 — > R. For the sake of simplicity, we will denote the induced function ||q|| H > B{q) 
by the same letter B, and define the radius R := supsuppB > 0. 

Due to the rotational symmetry of B, Noether's theorem implies the existence of a second 
constant of motion besides the kinetic energy. In order to find an explicit formula for it, 
we remark that (1) can be written as a Hamiltonian system in an alternative way. Since 
lu — luq = B{q)dq\ A dq2 = B{r)rdr A dip — da is exact with a = ( — J R B(r)r dr)dip 1 (f ) is 
generated by the Hamiltonian H(q,p) = l>\\p — A(q)\\ 2 with the magnetic vector potential 

A{q) := (y^ B(r)rdr ) j9 ' 

but now taken with respect to the standard symplectic form ljq. The corresponding La- 
grangian is given by 

L(q,q) = ±\\ q \\ 2 + (A(q),q). 

Proposition 2.1. The quantity I := (q + A(q), —Jq) is an integral of motion. In fact, the 
system (I) with a rotationally symmetric magnetic field B = B(\\q\\) is completely integrable. 

We call / the magnetic momentum. 

Proof. We apply Noether's theorem (see, e.g., Abraham and Marsden [AM]) to the one- 
parameter family 

h s : M 2 -j. M 2 , h s (q) - ( C ° SS q =: T s q. 

y ' \sins coss J 

For a solution q — q(t) we have jj^h s {q) = T s q and hence 

L(h s (q), f t h s (q)) = i||r s( ?|| 2 + (A(T s q),T s q). 

Since our choice of the magnetic potential A satisfies A(T s q) = T s A(q), we obtain 

L(h s (q),f t h s (q)) = l\\q\\ 2 + (A(q),q) = L(q,q). 

Now Noether's theorem yields that 

1 = (f (9,9), IjLo^)) = (q + A(q),-Jq) 

is an integral of motion. 

The energy is a second integral of motion. It is a straightforward calculation to show that 
both integrals are in involution, hence our 4-dimensional Hamiltonian system is completely 
integrable. □ 
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Now consider a solution q(t) of (1) entering and then leaving the support of the magnetic 
field B. At the unique entering and exiting points ? ± := q(f ± ) we have \\q ± \\ = R, hence 
A(q ± ) = 0. In view of Prop. 2.1, this yields the following result. 

Corollary 2.2. For an orbit entering and exiting supp-B, the entering and exiting angles 
arccos ((<j , (t ± ), Jq^) / (V2ER)) coincide. 

Note that in polar coordinates q — (r cos 9, r sin 9) the Hamiltonian satisfies H = \ (f 2 + 
r 2 9 2 ). Therefore we obtain, for some fixed energy E > 0, the expression 



I = r 2 t 



I B(r)r dr = ±ry / 2E - r 2 - / B(r)rdr, (2) 

J r J r 



showing that I does not depend on 9. 

In a constant magnetic field the trajectories are circles of a fixed radius, the Larmor 
radius. The question arises whether there are still circular orbits in a rotationally symmetric 
magnetic field. The first observation is that the curvature of a solution q of (1) at t equals 
^j^E^ ' Hence, the existence of a circular orbit in = H^ 1 {E) of radius r w.r.t. the 
origin of the configuration plane is equivalent to r satisfying the equation 

\B{r)\ 1 



V2E r 

This yields the following result. 

Lemma 2.3. For every energy E £ (Q,E°) with 

(B(r)r) 2 
E° := max - V 




r>0 2 

there are at least two circular orbits; for E — E° there is at least one circular orbit. 

For the following considerations we fix an energy E £ (0, E°\. If we denote by R + = 
R + (E) the largest radius of the circular orbits having energy E, then 

for r = R + 
for r > R + 

This outermost circular orbit plays an important role for the dynamics. Outside the disc 
of radius R + the magnetic field is too weak to capture orbits and prevent them from escaping 
to infinity. This is expressed by the following result. 

Proposition 2.4. Let E £ (0,E°] and xo = (qo,Po) £ with \\qo\\ > R + and (qo,po) > 0. 
Then there is a 5 > such that 

\\q(t)f> \Wo\\ 2 + St 2 

for all t > 0. 

Proof. On {||q|| > ||g ||} the function l^M ||g|j < 1 attains its maximum with maximal value, 
say, 1 — We consider the function f(t) := ^\\q(t)\\ 2 with derivatives 

f(t) = (q(t),p(t)) 

and 

f"(t) = 2E+(q(t),B(q(t))jp(t)) > 2E-\\q(t)\\\B(q(t))\V2E. 
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As long as ||g(f)|| > \\qo\\ the second derivative satisfies the inequality 

f"(t) >2E-V2E(l-^)V2E = 6. 

Since /'(0) > this holds for any t > and hence the claim follows. □ 

Remark 2.5. The calculation in the proof shows that t i— > ||g(£)|| 2 is convex while q(t) is 
outside the disc {\\q\\ < R + }. Hence, for an orbit staying outside the disc of radius R + , 
\\q(t)\\ either attains its minimum or converges to its infimum as t — > oo. 

Remark 2.6. The sign of B{R + ) ^ determines the orientation of the outermost circular 
orbit. The orbit winds clockwise around the origin if B{R + ) > 0, and counterclockwise if 
B{R+) < 0. 

For further investigations of the dynamics, especially of the orbits staying outside the 
i? + -disc, we make use of the integral and define /+ to be the value of / on the circular orbit 
of radius R + . For now we assume B{R + ) > and refer to Rem. 2.12 for the case B(R + ) < 0. 
We then have 

, rR 

I + = -R+V2E- B(r)rdr (3) 
Jr+ 

and call this quantity the critical magnetic momentum. 

In order to understand the motion of orbits entering the support of B, we consider the 
set of points 

UE = {(q,p)€V B \\\q\\=R,{q,p)<0} (4) 

through which orbits in g-space enter the ball of radius R. By taking the cosine of the angle 
between q and Jp we have Ue = S 1 X [—1,1]. Since (q, Jp) — r 2 9 we have a correspondence 
between 9 € [— "jp] and the entering direction in [—1,1]. Thus, orbits entering the 
support of B have, in view of (2), magnetic momentum / = R 2 8 e [—R\/2E, R\2E]. 

Lemma 2.7. If x = (q,p) G Ue has magnetic momentum I < I + then \\q(t,x)\\ > R + for 
all t e R. 

Proof. Outside suppi? the motion coincides with the free motion, hence the statement 
\\q(t, x)\\ > R > R + holds for all t < 0. Assume a trajectory intersects the circle of ra- 
dius R + . Then at the intersection point we have r = R + and, by uniqueness of the solution, 
r =/= 0. This implies 

[■R rR. 

I=(R + ) 2 9- B(r)r dr>-R+ ^2E - f 2 - B(r)rdr 
Jr+ Jr+ 

, r R 

>-R + V2E- B(r)rdr = I+, 

JR+ 

contradicting our assumption. □ 

Lemma 2.8. Let x = (q,p) 6 with \\q\\ > R + and I < I + . Then, for any p > \\q\\, there 
is a constant c p < such that 6 < c p as long as \\q(t,x)\\ < p. 
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Proof. The inequality 

• 1 



\q(t)\ 
1 



I(x) + [ B(r)rdr) 



P rP 



— Ti — TTITt' ( -R + V2E- / B(r)rdr + / B(r)rdr 



i 



fl+ ^||g(t)|| 



< 



k(*)ll 2 

R+V2E 



-R+V2E- / B{r)rdr 
Jr+ 



holds as long as \\q(t)\\ < p. □ 

Note that, in view of Lemma 2.7, points with I < I + cannot enter the disc of radius R + . 
Furthermore, for these trajectories we already know that assumes a global minimum or 

converges to its infimum for t — > oo. Which behaviour occurs can now be precisely described 
by the magnetic momentum. 

Proposition 2.9. For each (q,p) £ Ue we have 

1. 1 = 1+ ^ Mmt^Wq^W =R+. 

2. KI+ mint eR ||Q(t)|| > R+. 

Proof. As for the first assertion, we prove by contradiction. If does not converge 

to R + then it attains its minimum p > R + and, since Lemma 2.8 will assure 9 < 0, we have 

r-R 



I = -pV2E- [ B{r)rdr 
J P 

= -pV2E+ / B(r)rdr- / B{r)r dr 
Jr+ Jr+ 

, , /•« 

< -pV2E+ (p - R + )V2E- / B(r)rdr 

JR+ 

= -R+V2E- / B(r)rdr 

JR+ 



which contradicts the premise. The reversed implication is true since the curvature near 
{||g|| = R + } is negative and hence a trajectory with \\q(t)\\ — > R + has to have negative 
values of 8. Then the trajectory converges to the circular orbit, giving I = I + . 

For the second assertion, we have to show that cannot converge to its infimum 

as t — > oo. If this was the case we would have (q(t),p(t)) — > 0. By compactness there is a 
sequence t„ — > oo of times such that (q(t n ),p(t n )) — > (qoc,Poc) £ with (qoa,Poc) — 0. By 
Prop. 2.4 its trajectory leaves the support which contradicts the continuity of the flow. □ 

Example 2.10. We want to illustrate some of these results by an explicit example. Since 
the magnetic field B — B{r) need only be locally Lipschitz, we may consider the simplest type 
of Lipschitz continuous function having compact support, e.g., B : [0, oo) — > K where 



B(r) 



10(1 -r) if r£ [0,1] 
if r > 1 
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Then the magnetic momentum I for the energy E = 1/2 can be written as in (2). The graph 
of B, together with the graph of the function r 1— > 1/r, and f/ie fewe/ seis of I are shown in 
Fig. 2. 




(a) The graph of B (b) The level sets of I 



Figure 2: The magnetic field and its magnetic momentum 

Since there are only two intersections of the two graphs in Fig. 2a, we have precisely two 
circular orbits which are represented by the points (il^jO) in Fig. 2b; the outermost circular 
orbit has radius R + ~ .8873. The level set of the critical momentum is the level set where 
I — I + — I(R + , 0) ~ —.946. One clearly recognizes in this example, e.g., the characterization 
of orbits having magnetic momentum I > I + given in Prop. 2. 9. 

For points x £ Ue with I(x) < / + we now consider the angle 

Jo Jo h(v\r 

as well as, for I(x) < I + , the exit angle 

6 e (x) := 0(T(x),x) 

where T{x) denotes the exit time of x. Then Lemma 2.8 immediately yields the next result. 

Proposition 2.11. For fixed x the function 8(x, t) is strictly decreasing with respect to t. For 
x with I(x) = I + we have 6{x,t) — > — oo as t — > oo. We have e {x) — > — oo as I(x) — > I + . 

Remark 2.12. The previous observations hold for B{R + ) > 0. In fact, these calculations 
work in a similar way for the case B{R + ) < where one has to switch the sign of 9 since 
the circular orbit turns in the opposite direction. Thus, for B{R + ) < 0, we obtain 

, /•« 

I + = R + V2E- B(r)rdr (5) 
Jr+ 

as the critical magnetic momentum, and the analogue of Prop. 2. 9 now reads 

1. I = I + lim^oo ||g(t)|| = R + (no changes here). 

2. I> 1+ =^ mhi teR || 9 (t)|| >R+. 

Furthermore, we have 9 > for I > I + , and 9 e (x) — > oo as I{x) — > I + . 
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3 Symbolic Dynamics 



In the following we consider a magnetic field 

n 

B = Y J B k 

k=l 

where each component B k is rotationally symmetric with respect to its respective center 
q k £ K 2 , i.e., B k = B k {\\q - q k \\), with 

suppB fc C D k := {q £ R 2 | \\q- q k \\ < R k } and D k n A = for fc ^ Z. 

For a rotationally symmetric magnetic field B, we saw in Section 2 that, for energies below 
some threshold E° = E°{B) > 0, there is a circular orbit. In the case of several components 
we require the existence of circular orbits for each component B k and hence define 

E° := min E°(B k ) > 0. 

l<fc<n 

We point out that the magnetic momentum introduced in Prop. 2.1 is no longer a constant 
of motion since B is not rotationally symmetric anymore, but for each B k we may compute 
its local magnetic momentum I k (with q k taking the part of the origin) and obtain a local 
integral in the sense that I k is constant along a trajectory as long as it stays outside the 
other supports. 

For the following considerations we fix an energy E £ (0,E°] and let R+ = R k (E) 
denote the radius of the outermost circular orbit in D k . We write I k = lt(E) for the 
critical magnetic momentum for B k as in (3) or (5), depending on the sign of B k along the 
circular orbit of radius . 

Remark 3.1. Instead of looking at the value I{x) of the magnetic momentum we can consider 
the angle a between q and Jp. Since we are interested in certain entrance directions it 
is more convenient for the following to think of a critical angle instead of a critical 
momentum it . We have an orientation preserving homeomorphism between the momentum 
I £ [-RV2E, RV2E] and the angle a £ [-§ , §]. 

In what follows we will need two geometric conditions on the configuration of the compo- 
nents supp B k of the support of the magnetic field. Magnetic fields satisfying these conditions 
are said to be in general position, and we will derive a precise definition in the following (see 
Def. 3.2). 

1. The first condition will allow transitions from one support to any other. 

The easiest way to guarantee this would be to demand for the convex hull of two 
supports to have empty intersection with the other supports. 

However, this is a rather restrictive condition and we will use a weaker one, only 
demanding for certain parts of the convex hull — depending on the chosen energy — to 
have an empty intersection; see Fig. 3 for an illustration. 

2. The second condition will assure that we can choose an appropriate Poincare section 
which counts the revolutions around a given center q k in the right way. 

The transition corridors between two supports, as depicted in Fig. 3, depend on the 
orientation of their circular orbits, or equivalenty, on the sign of B along the orbit. For 
this we consider the four possible sign combinations. If for fixed k ^ I both signs are 
positive we consider the tangent line to dD k and dDi passing both supports to the left 
and the line which hits dD k with angle and dDi with angle see Fig. 4a. 
By A ky i and B\^ k we denote the part of the set of points on dD k and dDi between these 
two lines as shown in Figure 4b. If both signs are negative we basically consider the 
same lines as in the previous case but with left tangent replaced by right tangent. 
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Figure 3: Areas of empty intersection 




(a) a\ and (b) Ak,i and Bi : k 



Figure 4: Definition of and a[ , respectively, A^j and Bk,i 

If one sign is positive, say of Bk, and the other, say of Bi, is negative along the circular 
orbits we consider the left tangent to Dk which hits dDi with angle a z + and the line 
which hits dDk with angle and is tangent to D\ on the right hand side; see Fig. 5 
for a visualisation and the choice of Ak t i and Bi^ in this case. 

Let us point out that the sets Ak,i and B^k, as well as the critical angles and ou , 
depend on the energy value E G (0, E°]. We now give the precise definition of the geometric 
configuration we are going to impose. We denote by conv M the convex hull of a set M. 

Definition 3.2. Let B = X)fe=i Rk be a magnetic field on R 2 such the Bk are rotationally 
symmetric ( with respect to their respective centers qk ) and have pairwise disjoint supporting 
disks Dk = {q G R 2 | \\q — qk\\ < Rk}- We then say that B is in general position with 
respect to an energy E G (0, E°] if we have 

1. conv(j4fc,/ U 5/,fe) H supp B m = for distinct k, I, m G {1, ... , n} , and 

2- dD k \ {\J^ k A Kl U U/^fe ^ for all k. 

Remark 3.3. 1. Note that the Bk need not have a fixed sign. There may also be parts of 
the disc {\\q — qk\\ < Rk} where Bk vanishes. 

2. The first condition guarantees that all the transitions from Dk to any other Di are 
possible. The second condition allows us to put a Poincare section at an appropriate 
place; it will be a radial segment ending in the nonempty set considered there. 

3. Replacing the first condition by the much weaker condition conv(Z3/ £ U Di) n D m = 
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would still guarantee all the transitions, but exclude very simple configurations (n = 4 
is enough) which then would not satisfy Condition 2. 

Condition 2 assures the existence of points 

q* k £ 9£> & \(|JA M U|JS W ), 

l^k l^k 

and we pick such points ql and hold them fixed once and for all. With these points we now 
define the Poincare sections 

Pk,E = {{q,p) € £e I q = qk + \{ql - Ik) for some A £ (0, 1), (q - q k , Jp) ^ 0}, 
as well as the Poincare map 

n 

p: P E := |J P KE — > Pe"u {oo} 
fe=i 

by defining p(x) as the first point <&(t,x) € Pe for t > if such a point exists; we set 
p(x) = co if such a point does not exist. See Fig. 6) for illustration. Then p is continuous 
at x e Pe if p(x) G Pg. In particular, if x m — > e Pe and p(x rn ) j/oo € Pe, then the 
continuity of the flow yields p(xoo) = y^ and p is continuous at x^. 



p(x) = oo 

/ XL^^r ■ — i 


/ s Si * 







Figure 6: Definition of the Poincare map p 



Then 

A B = {x e P E I p*(a;) € Pe for all i e Z} (6) 
is an invariant set for p which, at this point, might still be empty. Let 

h: A E ^ S = {!,..., n} z (7) 
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denote the canonical coding map where h(x) = (si) such that p l (x) £ P Si for all i £ Z, and 
let er : S — > S be the left shift map, shifting a sequence (si) one position to the left. By 
construction we get the identity 

h o p ^ a o h. 

We now formulate the main theorem of this paper. 

Theorem 3.4. If the planar magnetic field B is in general position with respect to the energy 
E £ (0,E°], then there is a nonempty p-invariant subset A E C Ae such that h : A E — > S is 
continuous and surjective; in other words, p : A E — > A E is semi-conjugated to the full shift 
a:S ■ S. 

Since a dynamical system has at least the entropy of the one it is semi-conjugated to, 
we immediately obtain the following corollary for the corresponding magnetic flow. Here 
we use Bowen's definition of topological entropy (see, e.g., [Wa, § 7.2]), and denote the flow 
restricted to the energy surface £_g by <&£ := $|r x s e - 

Corollary 3.5. For n > 2 bumps and a planar magnetic field B = B^ in general 

position, there is a constant c > so that, for each E £ (0, E°], the topological entropy is 

bounded from below by 

hto P (<S>E) > cVE. 

Proof. By the monotonicity of the topological entropy with respect to the inclusion of in- 
variant sets [KH, Prop. 3.1.7(1)], it suffices to show that the topological entropy of the flow 
$_b e := $|rx_b e satisfies the same lower bound, with Be ■— 3>(R, Ae) and Ae C P E defined 
in (6). 

We now use Thm. 3.4 above. By the decrease of the topological entropy with respect to 
semi-conjugacies [Wa, Thm. 7.2], ht op (&B E ) is at least the topological entropy of the special 
flow ip given in [KH, §0.3] of the full shift with a certain roof function Te ■ S — » (0, oo). 
Indeed, as a surjection, the semi-conjugacy h from (7) has a left inverse S — > A' E , and 
choosing such an inverse, Te is defined as its pull-back of the Poincare return time 

Te ■ A' E — > (0, oo), f E (x) := inf{t > | $(t, x) £ P E }. 

We claim that there is a constant d > such that, for each E £ (0, E°), Te is uniformly 
bounded above by d / yE. This is true for points of A' E being mapped to a Poincare surface 
belonging to a different bump, since on the speed equals \/2E and the mutual distances 
of the configuration space projections of the Poincare surfaces are uniformly bounded from 
above. Then the Poincare time is bounded from above by that distance bound, divided by the 
speed. For points of A^ being mapped to the Poincare surface of the same bump, a similar 
argument applies, since the total length of the trajectory between successive intersections of 
Poincare surfaces is bounded from above by a constant times the circumference of the disk 
D k , 

By applying Abramov's formula, see Thm. 2.2.1 in [Si, Ch. 3], we finally conclude like in 
the proof of Thm. 14.4 of [Kn2] that 

d 

This implies the assertion of the corollary with c := I2IL™. □ 
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3.1 Proof of Theorem 3.4 



For the following considerations we need to define various sets and maps. Recall from (4) 
in Sect. 2 the definition of the set Ue of points in phase space through which g-space orbits 
can enter the support; due to the different centers of the supports this will now replaced by 

n 

U e = IJ U k ,E where U k ,E = {{q,p) € He | \\q - qk\\ = Rk, (q - qk,p) < 0}. 
fe=i 

We also need two additional maps that describe the relation between Pe and Ue- For some 
orbit entering the support let 

u : U E -> Pe U {oo} 

denote the first point $(t, x) G Pe, t > 0, to hit the Poincare section. For some orbit 
intersecting Pe let 

v : P E -> U E U {oo} 

denote the first point $>(t,x) € C/e, t > 0, where the orbit enters the support again. In case 
such points do not exist u(x), respectively v(x), are set to oo. 

In order to understand the basic mechanism to find a point that realises a prescribed 
itinerary, we start with a segment T C Uk consisting of phase space points entering Dk- We 
assume that the trajectory of one endpoint of F passes Dk tangentially to the left and the 
other endpoint has critical momentum i.e., its trajectory converges to the outermost 
circular orbit inside suppBfc, see Fig. 7a. 

Then we find in T a subsegment of points whose trajectories hit Pk exactly j times before 
leaving Dk towards Di. In particular, the resulting segment of points in Ui has the same 
configuration as the original segment entering Uk'- the trajectory of one endpoint passes Di 
tangentially to the left and the other endpoint has critical momentum This behaviour 
is illustrated in Fig. 7b and made precise in the following lemma. 




(a) Endpoints of original set (b) Endpoints of resulting set 



Figure 7: Shooting mechanism for j = 2 

Lemma 3.6. Let the magnetic field B be in general position with respect to the energy E. 
Let k 7^ / G {1, . . . , n} be different indices of supports and j G N. Let 7 : [a, b] Uk with 
h (7(a)) = -RkV2E and I k {f(b)) =/+.'" 

Then there exists [a, b] C [a, b] such that for all points x G im7|r- & in the image o/7|r s £1 
we have 

1. p l (u{x)) G Pk for < i < j — 1 

2. p»' («(»)) i p k 

3. v(pj{u(x))) G U 

and the curve 71 = (v o p> o u o -y)| ^ ^ : [5,6] — > Ui satisfies //(71(a)) = —Ri\j2E and 
7,(^(6)) =J+. 
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This description is valid for the case of the magnetic fields being positive along the 
outermost circular orbits, and the statement (as well as the proof) are also given for this 
case. The general situation is then treated in Rem. 3.7. 

Proof. Let q a , qb € dDk denote the two endpoints of dAk,i as shown in Fig. 8. 




Figure 8: Definition of q a and q b 



Since Jfc(7(s)) — > l£ for s — > b we have 9k(l(s)) — > —oo as s — > b by Lemma 2.11. Hence 
there are parameters s a < Sb in [a, b] such that the trajectory of 7(s a ) exits through q a , 
the trajectory of 7(s&) exits through qb, and for s a < s < Sb the forward trajectory hits Pk 
exactly j times before leaving Dk- The trajectory of 7(s a ) passes Di on the left hand side, 
and by Cor. 2.2 the trajectory of 'j(sb) passes to the right of the line hitting Di with angle 
af , hence it hits with some angle bigger than af , i.e., its magnetic momentum is bigger 
than 

Therefore there are new parameters a < b in [s Q ,Sf,]_ such that the trajectory of 7(a) 
passes tangentially to the left of D\, the trajectory of 7(6) hits Di with angle equal to , 
i.e., in a point with momentum equal to and for s£ [a, b] the trajectory of 7(5) hits Di. 

Then 71 : [a, b] — > Ui with 71 (s) = (v o p> o u o 7)(s) is well defined and possesses all the 
properties stated above. □ 

Remark 3.7. The mechanism in the proof of Lemma 3.6 is manufactured for the situation 
where both magnetic fields have postive values on their outermost circular orbits. We argue 
that an analogous result holds for all other cases as well. 

Indeed, in the case when the two magnetic fields have different signs on their outermost 
circular orbits, e.g., when > and P>i < 0, we can choose the curve 71 in such a way that 
the trajectory 0/7(0) passes tangentially to the right of Di, i.e., 7i(7i(a)) = Ri\/2E, and the 
trajectory 0/7(6) hits Di with angle , i.e., Ii(ji(b)) = 

Finally, an analogous mechanism works when both and B\ take negative values on 
their outermost circular orbits. 

The mechanism described above allows us to show that for every half-infinite sequence 
there is a point realising this prescribed itinerary. 

Proposition 3.8. For any half-infinite sequence (fc;) G {l,...,n} N ° there exists xq G Pk Q 
such thatp l (xo) G Pk { for each i>0. 

Proof. We start with a curve 7 : [0, 1] — > Uk such that /fc o (7(0)) = —Rk "v2E (respectively 
Rk V2E if Bk is negative along the orbit) and /fe (7(l)) = l£ . 

Let us assume first that (fc^) will not become constant eventually. Let j\ denote the 
number of consecutive values of fco, i.e., ki = kg for i < j\ — 1 and kj t 7^ ko- Then, 
by Lemma 3.6, there are a parameter interval [a^\b^] C [0,1] such that p l (u('j(s))) G 
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P ka = Pki for < i < ji — 1 and some corresponding curve 71 : [a^\ b^} — > U kji ■ Applying 
Lemma 3.6 to the curve 71 and the number j 2 of consecutive values of kj 1 , we get [a^ , b^} C 
[ a W,&W] such that p l (u(-/(s))) G P k]i = P ki for all s G [a( 2 \M 2 )] and ji < i < ji + j 2 - 1. 

Hence, by iteration, we get [d m \ M m )] C [a^"'- 1 ), fet" 1 " 1 )] such that p i (u('y(s))) G P fej 
for i (at least) up to m — 1. Then there is a point a; G f) m >o 7([a'" l \ fr'" 1 ^]) 7^ and hence 
Xq = u(x) G P ko satisfies p l (xo) G Pfc; for i > 0, which proves the claim if (fc^) does not 
become constant. 

If (hi) becomes constant eventually, we follow the same procedure but at the index where 
(ki) becomes constant, we simply choose the point with critical momentum and need not 
iterate further. □ 

In order to show that every bi-infinite sequence can be realised by some trajectory, we 
shift the sequence and show that the points obtained by Prop. 3.8 converge to a point which 
has the prescribed itinerary. Let us define 

A E := {x G A E I \\q(t,x)-q k \\ > R£ V t G R, k = l,...,n}} 

as the subset of points whose trajectories stay outside the open disks {\\q — qf-\\ < R k } for 
all times. 

In a first step, we need to show that h\ A > is surjective. Let (ki) G {1, . . . , n} z be given. 

For each m G N, Prop. 3.8 yields a point y m G Pk- m such that p l (y m ) G Pfc;_ m for i > 0, 
and we define 

x m := p m (y m ) G P ko . 

Since Pk is compact we have a convergent subsequence, which we denote by x m again, i.e., 
we have x m — > Xoo G P ko asm-> 00. We claim that Xoo lies in A E and satisfies h(xoo) = (ki), 
i.e., 

p 4 (xoo) e Pk> 

for all i G Z. 

In order to prove that Xoo G A B we first claim that the trajectory of a^oo does not intersect 
dP E . Arguing by contradiction, we assume that this happens at the index i = 0, say, so that 
we have 

Xoc = (goo, Poo) € dP k(l:E = {(q,p) € S B I q = q ko or q = q* ko or (q - q ka ,Jp) = 0}. 

Now, (qoc — qk i JPoc) — is not possible since (q — q k(n Jp) = \\q — qk a \\ 2 9 < by Lemma 2.8. 
This also implies that q = q ko cannot occur. Therefore, we must have qoo = g£ where three 
cases can occur, each of which will lead to a contradiction. 

Case 1: k\ 7^ k$. Then the g-space trajectories q(t,x m ) hit D kl , and by continuity 
of q and closedness of the support, so does q(t, Xoo)- But this contradicts the choice of 
fc ,fen as shown in Fig. 9. 

Case 2: fc_i 7^ k$. Similar to Case 1, the backward trajectory of Xoo in g-space would 
hit D k l , contradicting the choice of qfc ^ B ka ^ k l . 

Case 3: k-\ = ko = k\. Then \\q(t,x m ) — q ko \\ < R ko for \t\ < to with to > independent 
of m. Hence \\q(t, Xoo) - q ka \\ < R ka = || 9oo - g feo || for \t\ < t . Then ( 9oo - q kt)1 Poo) = 
2i^\\q(t, Xoo) — qk l| 2 | t=0 = 0, i-C-, the trajectory passes tangentially to D k() and satisfies 
\\q(t, Xoo) — qkoW > Rko for small \t\. But this contradicts \\q(t,x m ) — q ko \\ < R ko . 

Therefore the assumption Xoo G dP ko is false and, hence, the trajectory of Xoo does not 
intersect dP E at all. By continuity we finally get 

P l (xoo) = lim p l (x m ) G P ki 

m — > 00 

for all i G Z. This finishes the proof that h\ A > is surjective. 
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Figure 9: Contradiction in case 1 



It remains to show that h\A E is continuous. For this it is sufficient to show that the 
time between two consecutive intersections of the Poincare section is uniformly bounded. 
Here, two cases are possible, namely the intersections can occur in the same support, or in 
different supports. Lemma 2.8 gives a bound for the time how long a trajectory can stay 
inside the outer annulus in one support. Furthermore, outside the support the trajectories are 
straight lines and, since ||g|| = v2E, the time is proportional to the length of the trajectory. 
Therefore, we have a uniform bound for the length of the segments between two consecutive 
intersections of the orbit with the Poincare section and a uniform bound for the time. This 
proves the continuity of h on A E . 

Thus we have found a p-invariant set A E ^ on which h is continuous and surjective. 
This finishes the proof of our main result, Theorem 3.4. 

4 Conclusion and further remarks 

We have shown that a planar magnetic field with disjoint localized rotationally symmetric 
peaks creates positive topological entropy. To the best of our knowledge, this is the first 
instance where such systems have been investigated. In some sense, the classical potential 
analogue of our magnetic system should be the motion of a charged particle in a potential 
field with several Coulomb-type singularities. This system has been extensively studied in 
the literature, and it is well known that it carries positive topological entropy. 

We point out that the existence of a second integral, although we made some advantage 
of it in our proofs, is not really necessary. Indeed, we only need the rotational symmetry in 
the region outside the outermost circular orbit. Inside the discs of radii R~£, the magnetic 
field may behave as wildly as it wants; this will not affect any of our constructions. Actually, 
we expect that even outside the discs of radii we do not need any symmetry. The mere 
existence of a (hyperbolic) outermost closed orbit bounding a convex region in the q-planc 
should be enough to detect symbolic dynamics. This will be the topic of future research. 

Finally, the question arises whether there are higher-dimensional analogues of a magnetic 
field with localized peaks. At the moment, however, it is not clear what an appropriate 
generalization should be. Again, this will be investigated in the near future. 
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